Von Neumann Entropy Spectra and Entangled Excitations in Spin-Orbital Models 
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We consider the low-energy excitations of one-dimensional spin-orbital models which consist of spin waves, 
orbital waves, and joint spin-orbital excitations. Among the latter we identify strongly entangled spin-orbital 
bound states which appear as peaks in the von Neumann entropy (vNE) spectral function introduced in this 
work. The strong entanglement of bound states is manifested by a universal logarithmic scaling of the vNE with 
system size, while the vNE of other spin-orbital excitations saturates. We suggest that spin-orbital entanglement 
can be experimentally explored by the measurement of the dynamical spin-orbital correlations using resonant 
inelastic x-ray scattering, where strong spin-orbit coupling associated with the core hole plays a role. 
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Introduction. — The spin-orbital interplay is one of the im- 
ortant topics in the theory of strongly correlated electrons 
1J. In many cases, the intertwined spin-orbital interaction is 
decoupled by mean-field approximation, and the spin and or- 
bital dynamics are independent from each other. Thus a spin- 
only Heisenberg model can be derived by averaging over the 
orbital state, which successfully explains magnetism and opti- 
cal excitations in some materials, for instance in LaMnOs H. 
But in others, especially in t-2 9 systems yfl, the orbital degen- 
eracy plays an indispensable role in understanding the low- 
energy properties in the Mott insulators of transition metal 
oxides (TMOs), such as LaTi0 3 H, LaV0 3 and YV0 3 H, 
and also in recently discussed RbC<2 The well known 
cases are also strong spin-orbit coupling which leads to locally 
entangled states |7J|, and entanglement on the superexchange 
bonds in K3CU2F7 Jit]. For such models, the mean-field-type 
approximation and the decoupling of composite spin-orbital 
correlations fail and generate uncontrolled errors, even when 
the orbitals are polarized |Q. The strong spin-orbital fluctu- 
ations on the exchange bonds will induce the violation of the 
Goodenough-Kanamori rules 1 1 Oil - Furthermore, the flavors 
may form exotic composite spin-orbital excitations. 

Model and system. — A paradigmatic model derived for a 
TMO in Mott-insulating limit is the one-dimensional (ID) 
spin-orbital Hamiltonian, which reads 



H = -J 



Si ■ 5*1+1 +x) Tj • T i+ i + y 



(1) 



where S; and 2] are spin- 1/2 and pseudospin-1/2 operators 
representing the spin and orbital degrees of freedom located 
at site i, respectively, and we set below J = 1. It is proposed 
that ultracold fermions in zig-zag optical lattices can repro- 
duce an effective spin-orbital model Bill] . For general {x, y}, 
the model fl} has an SU(2)®SU(2) symmetry. An additional 
Z2 bisymmetry occurs by interchanging spin and orbital op- 
erators when x = y. In the case of x = y = \, Hamiltonian 
([TJ reduces to a SU(4) symmetric model, which is exactly sol- 
uble by the Bethe ansatz Il2lll3ll . There are three Goldstone 
modes corresponding to separate spin and orbital excitation, 



as well as composite spin-orbital excitations in case of J < 0, 
in contrast to a quadratic dependence of the energy upon the 
momentum in the long-wave limit for J > 0. The spectra 
of elementary excitations are commonly not analytically sol- 
uble away from the SU(4) point. We will, however, show 
that the low-energy excitations can be analytically obtained 
in some specific phases in the case when J > 0, and this 
offers a platform to study the spin-orbital entanglement. In 
this Letter, we go beyond the ideas developed for spin sys- 
tems n 1 411 . We demonstrate that spin-orbital entanglement en- 
tropy clearly distinguishes weakly correlated spin-orbital ex- 
citations from bound states and resonances by its magnitude 
and distinct scaling behavior. We propose how to connect the 
entanglement entropy with experimentally observable quanti- 
ties of recently developed spectroscopies. 

von Neumann entropy. — Currently, concepts from quantum 
information theory are being studied with the aim to explore 
many-body theory from another perspective and vice versa. 
A particularly fruitful direction is using quantum entangle- 
ment to shed light on exotic quantum phases E0. En- 
tanglement entropy even distinguishes phases in the absence 
of conventional order parameters II 171 . In general a many- 
body quantum system is subdivided into A and B parts and 
the entanglement entropy is the von Neumann entropy (vNE), 
<S v n = —Tr{ pa log 2 pa}, where pa = Tr b{p} is the reduced 
density matrix of the subspace A and p is the full density 
matrix. The vNE is bounded, 6> v n < log 2 dimp^, and easy 
to calculate. Experimental determination appears harder, yet 
there are proposals involving transport measurements in quan- 
tum point contacts iflitl . 

Interestingly the vNE scales proportionally to the boundary 
of the subregion obtained by the spatial partitioning Ill9ll . The 
dependence of the boundary or area law can be traced back to 
study of black hole phy sics lEoll and was extensively exploited 
for ID spin chains 112 ill . If the block A is of length I in a sys- 
tem of length L with periodic boundary condition, the vNE of 
gapped ground states is bounded as Si = 0(1), while a loga- 
rithmic scaling Si = clog 2 I + 0(1) (L > ! > 1) has been 
proven to be universal property of the gapless phases in criti- 
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FIG. 1 : (color online). Spin-orbital entanglement 5 v n in the ground- 
state of the spin-orbital model {I} as a function of x and y and system 
size L = 8. The (red) dashed lines mark the critical lines determined 
by the fidelity susceptibility (see text). The two-site configurations 
in phases I-IV are shown on the left. The two orbitals per site are 
degenerate (their splitting is only for clarity of presentation). 




FIG. 2: (color online), (a) Energy spectra of 40-site spin-orbital sys- 
tem at a; = y = 1/4. Dashed lines in the spin-orbital continuum 
denote the spin, orbital and OBS excitation, all degenerate; the (red) 
solid line below corresponds to the BS. (b) The decay rate V of the 
OBS for different momenta Q with y = x at L — > oo. 



cal systems by the underlying conformal field theory 112 211 . A 
violation of the area law is expected for the low-lying excited 
states of critical chains l23ll . To date, measurements of the 
vNE for subdivision of degrees of freedom other than in spa- 
tial segmentation have not been fully explored. In a composite 
system containing spin and orbital operators, the decomposi- 
tion of different flavors retains the real-space symmetries. 

Phase diagram. — A quantum phase transition (QPT) is 
identified as a point of nonanalyticity of the ground state and 
associated expectation values in the thermodynamic limit. To 
shed light on the phase boundaries, we first consider two sites 



|2J, H 12 = -\{S\ 2 - §1 -SI + 2x)(T? 2 - fl - T% + 2y), 



where S\ 2 = S\ + S 2 and T\ 2 = T\ + T 2 . A pair of spins 
(orbitals) can form either a singlet with S\ 2 = (T\ 2 = 0) or 
a triplet with Si 2 = 1 (T12 = 1), and various combinations 
of quantum numbers correspond to different phases shown in 
Fig. [TJ In phase I, the state with S12 = 1 = T\ 2 has the 
lowest energy, and thus the energy per bond is > e xy = 
— (x + l/4)(y + 1/4). For a larger system with L bonds, we 
have Eq(H) > Le xy . On the other hand, taking a ferro-ferro 
state |0) as a variational state, Eq(H) < Le xy . Therefore, 
the energy of phase I is exactly Eq(H) = Le xy and the ferro- 
ferro state is the corresponding ground state. 

Without prior knowledge of order parameter, various char- 
acterizations from the perspective of quantum information 
theory can be used to identify phase boundaries. One often 
used tool is the vNE 12511 . Tracing orbital degrees of freedom, 
we obtained the spin-orbital vNE <S v n for the ground state of 
L = 8 chain in the Hilbert subspace of S z = T z = 125 1. 
However, here we find that the vNE of the ground state does 
not distinguish phase I from phase II or IV — all three phases 
having 6> v n = (see Fig. [T). Therefore we use the quantum fi- 
delity to quantify the phase diagram [26J] . The fidelity defined 



as follows, T(X,8X) = |(* (A)|*o(A + SX))\, is taken along 
a certain path {x(X), y(X)} and reveals all phase boundaries. 
The fidelity susceptibility, xf = ~ (2 In J 7 ) / (5X) 2 \s\^o, ex- 
hibits a peak at the critical point, and can be treated as a ver- 
satile order parameter in distinguishing ground states 112711 . It 
signals the phase boundaries shown in Fig. Q] Remarkably, 
the phase diagram found from the fidelity susceptibility for 
larger systems is the same as the one for L = 2. 

Excitations. — In phase I of Fig. Q] with boundaries given 



by: x + y 



and y = —\, the spins and orbitals 



are fully polarized, and the ferro-ferro ground state |0) is dis- 
entangled, i.e., can be factorized into spin and orbital sector. 
It is now interesting to ask whether: (i) the vanishing spin- 
orbital entanglement in the ground state implies a suppression 
of joint spin-orbital quantum fluctuations, and (ii) collective 
spin-orbital excitations can form. Using equation of motion 
method one finds spin (magnon) excitations with dispersion 
w s(q) = (1 + y) (1 — cos q), and orbital (orbiton) excitations, 



Ut{q) = (\ + x)(l — cos q) 12811 . The stability of the orbitons 
(magnons) implies that x > — | (y > — 7), and determines 
the QPT between phases I and II (IV), respectively, while the 
spin-orbital coupling only renormalizes the spectra. 

For our purpose, it is straightforward to consider the prop- 
agation of a pair of magnon and orbiton along the ferro-ferro 
chain, by simultaneously exciting a single spin and a single 
orbital. The translation symmetry imposes that total momen- 
tum Q = 2rmr/L (to = 0, • • ■ , L — 1) is conserved during 
scattering. The scattering of magnon and orbiton with initial 
(final) momenta - q, § + q} + <?'}) and total 

momentum Q is represented by the Green's function i29ll . 
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for a combined spin (S Q ) and orbital (T Q ) excitation. 
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FIG. 3: (color online). The vNE distribution of 400-site spin-orbital 
system in subspace Pst = 1 for x = y = 1/2 and different momenta 
Q. Isolated vertical lines indicate the BS, with dispersion given by 
the (red) dashed line. The OBS in the center of spectra is damped. 
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The spin-orbital continuum is given by 0(Q, q) = 
q) + u>t{® + q). In the noninteracting case, the Green's 
function exhibits square-root singularities at the edges of the 
continuum l3(ill . Due to residual, attractive interactions spin- 
orbital bound states (BSs) are shifted outside the continuum 
l32ll . see Fig. EJa). The collective mode is determined 
i Mcos % - cos qfl [u - fi(Q, q)] = 0. The 



I24j,|31l 

by 1 + 

analytic solution of this equation is tedious but straightfor- 
ward. The collective BS with dispersion ujbs(Q) is well sep- 
arated from the spinon-orbiton continuum [Fig. |2ja)] at large 
Q. In the long-wave limit the BS energy coincides with the 
Arovas-Auerbach line 0311 . i.e., the boundary of the contin- 
uum, yet the BS remains undamped for x + y > \. 

In addition, a collective mode of spin-orbital resonances, 



l*(Q)> 



m+Zl 



(3) 



m,l 



occurs inside the continuum. Here < I < L — 1 de- 
notes the distance between spin and orbital flips. Remark- 
ably, the spin and orbital flips are glued together at the same 
site with ai(Q) = S^o at the SU(4) point [28]. This cou- 
pled on-site BS (OBS) is a coherent superposition of local 
modes, all of them with equal weight. It has dispersion 
^obs(Q) = x + y — i cosQ, which is degenerate with both 
uj s (Q) and u>t(Q) at x = y = 4, see Fig. |2ja). This is remi- 
niscent of the degeneracy of the three Goldstone modes at the 
SU(4) point for J = — 1 11121 1 1 311 - Moving away from the 
SU(4) point, the OBS decays due to residual interactions into 
magnon-orbiton pairs, and the mean separation £ of spin and 
orbital excitations increases, i.e., ai(Q) ~ exp(— //£), lead- 
ing in the thermodynamic limit to a finite linewidth defined 
by T = ImG -1 (Q, w) 13411 . The decay rate of the spin-orbital 



OBS increases with growing x > \ and also for decreasing 
momenta Q, as seen in Fig. Gfb). 

Entropy spectral function. — To investigate the degree of 
entanglement of excited states, we introduce the vNE spectral 
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FIG. 4: (color online), (a) Scaling behavior of entanglement entropy 
5vN of the spin-orbital BSs for Q = 0.87T. Lines represent loga- 



rithmic fits <Svl 



log, L + cq, with Co : 



-0.659, -1.059, -1.251, 



respectively, (b) The scaling behavior of entanglement entropy of the 
OBS for x = y — 1/2. Lines are fitted by 5 v n = c\/L + Co, with 
c (ci) = 3.69 (380.5), 3.37 (138.4) and 3.31 (47.6) for Q = 0.8tt, 
0.6-7T and 0.2-7T. The inset shows the logarithmic behavior of <S v n for 
the OBS with Q = 0.8rr and x = y = 1/4. 



function in the Lehmann representation, 

Svn(Q, w) = - Trip^ 5 log 2 P M}6 {lo - oj n (Q)} , (4) 

a 

where (/x) = (Q,u> n ) denote momentum and excitation en- 
ergy, and = Tr |5 , „(Q))(\l/ n ((5)| is obtained by tracing 
the orbital degrees of freedom. Let us first consider the sym- 
metric case, i.e., x = y. The Hilbert space can be divided 
into two subspaces characterized by the parity Pst of the in- 
terchange of S <-> T, which is odd or even. Translation sym- 
metry allows us to express the reduced density matrix p s in a 
block-diagonal form, where each block corresponds to an irre- 
ducible representation labeled by total momentum Q and par- 
ity of exchange symmetry Pst- The vNE can be obtained by 
diagonalizing separately these blocks. In particular, the non- 
degenerate eigenstates with odd parity can be explicitly cast 
in the form 4= (SZ ,„ Tr,„ . — 5Z,„. TZ, n )|0). Con- 

V2 K Q/2-9 Q/2+9 Q/2+9 Q/2-qJl ' 

sequently, the singlet-like pair results in <S v n = 1. For other 



spin-orbital eigenstates with Pst = 1, <Sl > 1, except the 
pure spin and orbital waves. Interestingly, we find that the 
parity is still conserved in subspace Q = for x ^ y. The 
strongly entangled spin-orbital BSs are reflected by peaks in 
the von Neumann spectra 5 v n(w), shown in Fig. [3] As mo- 
mentum Q decreases, the OBS-peak in the center of spectra 
gets broader, implying a shorter lifetime. 

Inspection of vNE spectra shows that the entanglement 
reaches a local maximum at the BSs. Finite size scaling 
of vNE of spin-orbital BSs reveals the asymptotic logarith- 
mic scaling £> v n = l°g2 L + cq shown in Fig. |4|a). The 
same logarithmic scaling is found for the OBS at the SU(4) 
point x = y = j, as seen in the inset of Fig. HJb). How- 
ever, far away from the SU(4) point the scaling is entirely 
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FIG. 5: (color online). The spectral function of the on-site excitation 
Ao{Q,uj) for: (a) x = y = 1/4, (b) x = y = 1/2; (c) the nearest- 
neighbor Ai+(Q,w) for x = y = 1/2. The momenta range from 
7r/10 (bottom) to 97r/10 (top); the peak broadening is r\ = 0.01. 
Dashed (red) and dotted (green) lines correspond to the BS and OBS, 
while gray dash-dot lines indicate the boundaries of the continuum. 



different and the entropy of the OBS scales as a power law, 
SvN = ci/L + cq, as seen in Fig. SJb). This change of scaling 
from logarithmic to power law in 1/L is controlled by the cor- 
relation length £ measuring the average distance of spin and 
orbital excitations in the OBS wave function OJ. From Eq. 
© and ai(Q) ~ exp(— //£) we obtain, 



S vN ~log 2 {L/(l+£)}, 



(5) 



which yields log 2 L at x = y = 1/4 where £ = 0. As £ 
increases the correction to the vNE is oc — log 2 (l + £). Far 
away from the SU(4) point, the OBS is damped and £ becomes 
extensive, i.e., £/Z w Co — ci/Z, and the vNE approaches a 
finite value with a correction oc 1/L as shown in Fig. 0b). 
This close correspondence of the vNE of bound states and the 
correlation length £ suggests to use the dynamic spin-orbital 
correlation function as a probe of spin-orbital entanglement 
and as a qualitative measure of the vNE spectra. 

Spectral functions. — Returning to TMOs, one realizes that 
joint spin-orbital excitations are not created in the ferro-ferro 
ground state in photoemission spectroscopy because of spin- 
conservation. On the contrary, resonant inelastic x-ray scatter- 
ing (RIXS) |35|-[39l] is in principle able to measure the spectral 
function of the coupled spin-orbital excitations at distance I, 



— lim Im 

7T V->0 



,(0t 



1 



■>(0 



Eo-H-iri 



Here V 



(o) 
Q 



7i£ 



3 °3 ± 3 



(6) 

is the local excitation op- 



erator for an on-site spin-orbital excitation. We use as well 



,(i±) 



-i-yv 



JQj 



(S- 



1 ± s j 



-1)27 



for the nearest- 



neighbor excitation. In the RIXS process an electron with spin 
up is excited by the incoming x-rays from a deep-lying core 
level into the valence shell. For the time of its existence the 



core hole generates a Coulomb potential and a strong spin- 
orbit coupling that allows for the non-conservation of spin. 
Next the hole is filled by an electron from the occupied va- 
lence band under the emission of an x-ray. This RIXS pro- 
cess creates a joint spin-valence excitation with momentum 
Qin — Qout and energy u- in — w out , which can unveil the spec- 
tral function of the spin-orbital excitation. 

The on-site spectral function Aq(Q,uj) shown in Figs. 
Ha,b) highlights the OBS. At the SU(4) point [Fig. Ha)] it ap- 
pears as a (5-function, Aq(Q, uj) = 5{uj — wobs(Q)}, whereas 
in Fig. |3b) the OBS is damped and its intensity decreases 
strongly with Q. In the latter figure the BS at the low en- 
ergy side of the continuum appears as weak additional feature, 
while it is absent in (a), i.e., at the SU(4) point. The nearest 
neighbor spectral function Ai+ (Q, u>) in Fig. [5jc) shows both 
the spin-orbital continuum and the BS outside of the contin- 
uum. Notably, comparing with the vNE spectral function in 
Fig. [3] we find the same characteristic energies and similar 
intensity features as in the RIXS spectra. The spectral func- 
tion provides information of various correlations, which are 
ingredients to derive the reduced density matrices [40]. 

Summary. — In this Letter, we study a spin-orbital system 
and extend the analysis of entanglement to excited states by 
introducing the vNE spectral function. Our study demon- 
strates that even in cases where the ground state of a spin- 
orbital chain is fully disentangled, e.g., in the ferro-ferro state, 
(i) the spin-orbital excitations are in general entangled, (ii) 
maximal spin-orbital entanglement occurs for BSs which ap- 
pear as sharp peaks in the vNE spectra, and (iii) the vNE of un- 
damped BSs exhibits a logarithmic dependence on the chain 
length L. We propose to study the dynamic spin-orbital cor- 
relation function as a qualitative measure of the vNE spectra, 
and suggest to use here RIXS as a promising technique. 
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